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NAVAL  ENGINEER. 

A  graphical  method  of  determining  normalized  curvature 
patterns,   deflection  patterns,   and  frequencies  for  the  natural 
modes  of  interest  of  a  ship  treated  as  a  non-uniform,   free-free, 
transversely  vibrating  beam  is  presented.     Corrections  for  shear 
deflection  and  rotatory  inertia  are  integrated  into  the  method, 
which  is  based    on  procedures  developed  by  Stodola  and 
Schadlofsky.     Integrations  are  performed  mechanically  by  an 
integraph.     The  method  is  demonstrated  by  determining  the  first 
three  horizontal  modes  of  a  three -eighths  scale  model  of  a 
submerged  submarine. 

The  method  is  found  to  be  convergent,    yielding  acceptable 
results,   and  taking  account  of  the  secondary  effects  considered. 
Three  cycles  of  iteration  for  each  of  the  first  three  non-trivial 
natural  modes  of  the  model  yielded  frequencies  of  417,    838,   and 
1391  cycles  per  minute,    respectively.     Patterns  obtained  are  tabu- 
lated in  detail. 

The  process  developed  is  recommended  for  consideration  by 
small  research  facilities  with  limited  budgets. 
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NOMENCLATURE 

A       =    total  cross -section  area  of  hull  structure  considered  effective 
in  betiding  and  shear,    (sq.  inche) 

*    axis  translation  of  assumed  or  calculated  deflection  for  ortho- 
gonality with  rigid  body  translation 

B       *    rotation  coefficient  tor  axis  of  assumed  or  calculated  deflection 
lor  orthogonality  with  rigid     ody  rotation 

C,  C    . .«  adjustment  coefficients  Lor  correcting  assumed  or  calculated 
deflections  for  orthogonality  with  previous  modes 

G.      «  constant  of  integration  applied  to  shear  (equal  to  aero) 

m 

CL      *  constant  of  integration  applied  to  bending  (equal  to  aero) 
C        *  constant  of  integration  applied  to  slope 

C.      «   constant  of  integration  applied  to  deflection 

6 

£        a   Young's  modulus,    I  K 

G        *    shear  modulus,     11. 5    x    10      psi 

I,  I-_      *  moment  of  inertia  of  hull  structure  about  vertical  axis  through 
centroid  of  each  section. 

K       ■  ratio  of  average  shear  stress  to  maximum  shear  stress  over 
a  cross -section  for  horizontal  transverse  shear. 

L       -  length  overall  of  ship,   ft. 

M      a  bending  moment,  ■%•!•« 


N       «  an  operational  constant,     »    /- — - —   m £  S* 
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KQMENCLVIUBE  (Continued) 
operational  constant,    *     J  J    i-*-~ —     m/  <*f  "T 


dx 

L        L 


O       «  shear  force,  lb 

T       •  period  of  vibration*   milliseconds 

W      *  an  operational  constant    «       jY,1'  m$  dx/L 

■  an  operational  constant    *      J  yf  mo    dx/L 


T  •  Y _.*       deflection  patterns  for  a  uniform  prismatic  bar  obtained  from 
consideration  of  simple  bending  theory  only,  according  to  Lord 
Kayleigh. 

f        ■       frequency  of  vibration,  cps 

g       *       force  of  gravity,   32.2  ft/sec 

m      *       virtual  mass  ptir  unit  length,  including  structure,   ballast,   fee- 

flooding  water,  etc. ,  plus  water  entrained  in  horizontal  vibration 

(lb  wTflr) 

w       *       mo  ,   virtual  weight  per  unit  length  (lb/ft) 

x       ■       coordinate  parallel  to  longitudinal  axis  of  ship,   measured  from 

bow  (ft) 
x       »       coordinate  of  longitudinal  centroid  of  virtual  mass 

y       a       coordinate  parallel  to  beam  of  ship,   measured  from  vertical 

plane  through  centroid  of  ship;    measure  of  horizontal  deflection 

y        •       the  deflection  pattern  obtained  from  a  cycle,   uefore  orthogona~ 

ligation 
y*     ■       value  of  unnormalizcd  deflection  curve  taken  at  the  normalizing 

point 
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NOMENCLATURE  Continued 

Z       »       coordinate  parallel  to  vertical  axis  of  ship  measured  from 
plane  through  centroid  of  ship. 

4        »     virtual  displacement  (pounds);  see 

(9  *  component  of  slope  due  to  shear  deflection 

T  *  component  of  slope  due  to  bending 

9"  *  component  of  slope  due  to  bending  before  axis  determinate*)  T^C^mf 

>1  *  deflection  before  determination  of  axis;  ji   +    c^   «*     y 

*J  «  slope  before  determination  of  axis;   r?  -f     c     «   y* 

jf  *  dummy  variable,   replacing    x 

6>  *  circular  frequency,    (radians/ second) 


Note:    Subscript  to  deflection  indicates  mode  number. 

Hyphenated  subscript  indicates  mode  number  and  Iterative  cycle  num- 
ber. 
Prime  indicates  differentiation  with  respect  to  x. 
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I.        INTRODUCTION 
.     Historical 


For  th*  last  seventy  years,   the  subject  of  vibrations  of  ship  hulls 
has  attracted  ever  increasing  attention  from  the  naval  architect.     Lariy 
investigators,  notably  Otto  Schlick,  emphasised  the  empirical  treatment; 
out  so  complex  is  the  Structure  of  a  large  ship,  and  so  varied  its  design 
features,   that  it   >ecam«  clear  that  this  approach  would  not  lully  B9rv  the 
need.     Ultimately  then,  we  must  look  to  the  theoretical  approach  for  the 
development  o:  adequate  formulations. 

Basically,   the  ship  is  considered  a  slender  free-free  elastic  beam, 
vibrating  transversely  to  its  long  axis,  for  which  the  theory  stems  from 
Lord  Mayleigh.    [©T*  * "rom  energy  considerations  he  developed  the  us  sic 
equation  based  on  simple   .ending.     He  then  showed  that  for  various 
assumed  mode  shapes,   the  one  for  which  the  energy  balance  gave  the  lowest 
frequency  was  the  correct  shape.     He  also  published  a  modified  equation, 
accounting  ior  the  effect  of  rotatory  inertia,   but,   considering  this  e;  ect 
negligible,  did  not  emphasise  It, 

noshenko  formulated  the  shear  flexibility  term  for  *  lyleigh's 
basic  equation,  and  thisi  set  down  the  well  known  limoshenko  mechanism, 
which  considers  the  contribution  of  bending,    shear,   and  rotatory  inertia. 
While  thi*  formulation  is  the  most  complete  available,   it  is  complex  and 
difficult  to  solve  exactly,   even  for  a  uniform  oeam,   and  in  recent  years  a 
numuer  of  papers  have  oeen  devoted  to  this  solution,  generally  involving 
Laplace  transforms  or  series  .-solutions.        11,    33,    34,   63,   et. al. 

Interest  in  the  non-uniform  u«&m  inspired  work  along  different  iines, 
since  a  beam  with  a  general  distribution  of  mas?,   inertia,  and  section  area 
is  not  amenable  to  exact  »oiution.     oue  approach  is  to  attem:  ate 

one  or  more  of  these  distributions  into  expressions  suitable  for  analytic 

*    Numbers  in  brackets  refer  to  references  listed  in  Appendix  F. 
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treatment  by  operational  calculus,   e.g.  C?»   10,    33,    37]   .     However,   it 
is  evident  that  a  reasonably  accurate  method  ox  solution  which  is  unre- 
stricted in  its  application,   i.e.,   can  tre*t  completely  general  distribution* 
of  structural  characteristics,  has  advantages  over  any  isolated  analytical 
solution,   however  rigorous.     Thus  we  shall  consider  the  iterative  methods, 
both  graphical  and  tabular. 

Stodoia's  work  on  turbines    [8}  suggested  treating  non-uniform  beams 
integrating  i  ayleigh**  differential  equation  numerically  or  graphically. 
His  methods  also  included  the  use  oi  the  funicular  polygon,  and  ■ a. yi sigh's 
rule  oi  minimum  frequency.     In  using  the  successive  integration  method, 
he  observed  that  while  the  two-noded  mode  behaved  well,   a  three -noded 
starting  assumption  diverged  during  iteration  until  it  returned  to  the  two- 
noded  mode,    he  failed  to  recognize  the  mutual  orthogonality  of  normal 
mode!!..     Viaueiio  introduced  this  method  into  general  use  in  the  field  of 
engineering,  and  J.    ..    iayior  and  Schadlofeky  applied  it  to  ships  (illustrated 
in  £6fc}  ,   as  is  lewis's  tabular  method).     All  worked  primarily  with  the 
first,  or  two-noded  mode,   and  employed  the  simple  -sending  iorm  of  the 
equation  oi  ^ayleigh,   though  lay  lor  recognised  the  importance  8  •* 

shenko's  shear  term  in  ship  vibrations.    Koch  [2]  proved  the  convergence 
oi  the  method,  and  more  important,  demonstrated  the  necessity  oi 
rigorously  maintaining  orthogonality  with  all  previous  modes  throughout 
the  iterative  method. 

Murray  demonstrated  the  importance  of  transverse  shear  in  ahip 
vibration*  £7u],   and  introduced  this  effect  into  achadloisky's  method, 
(rather  than  applying  it  as  a  subsequent  correction  to  frequency  as  did 
McQoldrick  £6*0  ).     This  refined  method,   first  employed  by  is   I         ,   and 
the  authors  in  work  on  a  submerged  submarine  [71]  ,   was  uoted  also  tor 
a  method  of  determining  the  slope  axis  by  integration  rts  of  the 
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expression  for  orthogonality  with  rigid  body  rotation  (pitch  and  yaw),   and 
lor  calculation  of  higher  mode*  than  the  first*   facilitated  by  application  of 
Koch's  intermodc  orthogonality  conditions. 

Other  graphical  methods  have  been  employed  oy  Southwell,   Gurnbel, 
Morrow,   and  other-.,    Uut  were  suitable  generally  for  the  first  mode  only. 

Beceotiy  the  principles  employed  in  these  iterative  methods  have  oeen 
applied  to  lump  parameter  formulations  for  solution  Ly  sequence  controlled 
digital  calculators  or  t>y  electrical  analogy.  £56,   67,   69]  .    In  particular, 

•  hi  -ind  Myklestad  developed  finite  difference  equations  for  use  in  the 
I  ay  lor  Model  Basin  digital  computer,   taking  account  of  simple  bending, 
shear  deflection,   and  rotatory  inertia.     Calculations  made  b,    this  method 
for  vibration  in  the  vertical  plane  for  the  attack  transport  U.S.  S.   NIAGARA 
with  and  without  che  secondary  effects  (shear  flexibility  and  rotatory  inertia) 
were  compared  with  experimental  data.     I  he  increasing  importance  of 
secondary  effects  with  increasing  frequency  was  demonstrated,   but  toff  this 
vessel  rotatory  inertia  was  relatively  unimportant  in  the  vertical  plane. 

B.      ..    active 

i  his  investigation  is  concerned  specifically  with  transverse  vibration 
in  the  horizontal  plane  of  a  submerged  submarine,     hecent  interest  in  the 
oodily  response  of  a  submarine  to  dynamic  loading  such  as  that  accompanying 
a  non-contact  underwater  explosion  has    ltd  to  the  development  of  methods 
of  modal  analysis.     These  methods,   and  other  demands  oi  modem  shi? 
design,   have  underscored  the  inadequacy  of  older  established  methods  of 
determining  natural  .  .odes  of  response  of  a  ship. 

In  particular,   then,   the  object  of  this  investigation  was  to  establish 
and  demonstrate  a  method  of  determining,   as  accurately  as  possible,   the 
curvature  pattern,   deflection  pattern,   and  natural  frequencies  oi  the  normal 
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.nodes  of  interest  la  a  transversely  vibrating  submarine  of  general  mass, 
inertia,   and  section  area  distribution.     Further,   the      method  was  to  uc 
capable  of  execution  without  resort  to  such  expensive  mechanisms  as 
digital  computers. 

C.    Method: 


Ihe  method  arrived  at  is  basically  similar  to  that  of  Schadlofsky 
L*b*Qi   involving  iterative  cycles  operating  on  an  assumed  deflection 
pattern,   each  cycle  comprising  the  four  successive  integrations  suggested 
by  iiayleigh's  equation.     The  Timoshenko  shear  deflection  correction  is 

ade  to  the  slope  o   caicd  from  bending,   and  the  rotatory  inertia  correc- 
tion of  Mayleigh  is  made  to  the  derivative  of  the  moment.     Koch's  ortho- 
gonality checks  are  made  to  filter  out  components  of  all  lower  modes,   and 
Murray's  t&sjpa  axis  determination  is  employed.     Integration  is  performed 
by  integraph  (Plate  I),   and  calculations  oy  commercial  desk  calculator. 
The  method  is  demonstrated  for  the  first  three  normal  modes,   using 
structural  data  from  an  existing  three -eighths  scale  model  of  &  submarine, 
including  an  allowance  for  entrained  water. 

i-hear  lag,   and  the  possible  inter&ctiou  between  shear  deflection 
and  rotatory  inertia  are  not  considered. 

It  is  noted  t  at  the  problem  is  formulated  as  an  eigenvalue  prob- 
lem,  in  widen  solutions  exist  only  for  a  discrete  set  of  squared  frequencies, 
each  of  which  corresponds  to  a  particular  mode  pattern,   except  that  the 
zero  eigenvalue  indicates  both  translational  and  rotational  vibration  of 
the  rigid  sody. 
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II    P  K  o  C  £  D  b  .t-  JL 

1  he  basic  theory  of  vibration  of  a  non-uniform,    transversely 
vibrating  beam  as  formulated  by  L,ord  Wu  yleigh,   considering  only  the 
simple  beading  effect,   can  he  expressed  briefly  by: 


d*  (  U^)     *      .    -nyco*  (i) 


where  the  assumption  oi  simple  harmonic  motion  ot  circular  frequency 
has  been  incorporated.     Four  successive  integrations  with  respect  to 
thin  equation  comprise  the  framework  upon  which  the  iterative  method  oi 
this  investigation  hat  developed.     The  deflection  and  irequency 

included  in  the  inertia  loading  term  on  the  right  side  ox  (1)  arc  assumed 
values,  and  the  deflection  indicated  on  the  left,   obtained  by  the    quadruple 
integration,   ie  the  first  estimate  oi  the  true  deflection  pattern.     Ihe  first 
estimate  of  the  square  of  the  true  circular  irequency  is  obtained  uy 
comparison  of  the  assumed  and  derived  deflection  antipodal  amplitudes,   i.e., 

<   •**  »liMl    »_j,  fj*)    »urt  (I) 

(y*)    final 

Note  also  that  the  second  integral  of  inertia  loading,  uMch  is  the  curve  of 

bending  moment,    MMl  Km  divided  by  the  curve  o.  local  ilexurni    ri 

before  the  last  two  integrations  are  performed,    since  Li  is  a  variable  in  x. 

We  first  modify  the  process*  by  applying  to  ihe  shear  curve  yielded  by 
the  first  integration  a  correction  for  rotatory  inertia,   as  follows: 

fUt  «        Q       -  JmT  (3) 

dx  A 

where  f  is  the  bending  component  oi  slope.     Since  this  quantity  occurs 

later  in  any  iterative  cycle  than  M,   we  must  use  the  Y  value  from  the 

previous  cycle  as  a  best  approximation,     ihus  this  correction    is.  not  made 
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until  the  second  cycle  of  each  mode. 

1  he  second  modification  of  the  process  is  the  correction  of  the  slope 
carve  obtained  b    integrating  the  curvature  distribution.     This  correction 
introduces  the  effect  of  shear  flexibility,   and  the  total  slope  resulting  is 
expressed  by: 

rfnr+fi  (4) 


where 


KG 


Here  the  shear  value  Q   has  been  obtained  before  the  slope  calculation, 
and  the  correction  can  be  made  in  each  cycle. 

Application  of  boundary  conditions  and  orthogonality  conditions, 
including  orthogonality  (with  respect  to  mass)  with  the  ^ero  eigenvalue 
motions  and  with  hll  modes  below  the  one  being  obt&ined,  at  the 
appropriate  points  in  and  after  each  cycle,   determine  the  four  constants 
of  integration,  and  refine  the  pattern  obtained  from  the  iterative  cycle. 

In  this  love*  tigs  lion,    since  mass  and  flcxurai  and  shear  rigidities 
have  general,   hence  irregular  distributions,   analytical  integration  was 
not  possible.     I  he  process  was  therefore  carried  out  in  each  instance  by 
most*   of  a  mechanical  integraph  with  a  Z-i"  carriage  beam  (.:.«mufuctured 
in  Switzerland  fej  G.   Coredij.     Ibis  instrument  draws  an  integral  curve 
when  its  stylus  is  made  to  trace  the  argument  curve  (Hate  I).     All  calcu- 
lations have  been  carried  out  and  plotted  at  intervals  of  one-fortieth  of 
hip  length,  as  this  was  found  to  be  the  greatest  interval  which  would 
detect  important  variationa  in  the  functions  used. 

The  sequence  of  operations  is  as  follows: 
(1)     Ihe  structural  characteristics  needed  are  computed  from 
available  plans  and  tabulated  (Appendix  E), 
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(2)  A  tabulation  of  functions  of  structural  characteristics  and  ot 
the  ccntroid  of  virtual  mass  which  will  be  used  throughout 
the  work  is  prepared.     (Table  IV,   Appendix  C).     lo  this 
table  are  added  a  tabulation    of   a  function  of  each  modal 
deflection  pattern  obtained,  for  use  in  orthogonali;eiag  tl 
sequent  modes. 

(3)  A  tabulation  of  constants  computed  for  repeated  use  U  prepared 
(Table  V,    Appendix  C). 

(4)  A  tabulation  of  the  &ayleigh  mode  patterns  (obtained  from 
application  of  his  simple -oending  theory  lor  a  uniform  slender 
beam)  is  prepared,   for  all  mode*  lo  te  obtained.     These  are 
normalized  patterns. 

(5)  The  .^.ayicigh  pattern  for  the  mode   -eing  calculated  is  orthogona- 
J.ized  with  translation,    rotation,  and  all  previous  modes  (there 
would  oe  none  for  the  two-aoded  mode).     In  practice,   the  Rayieigh 
pattern  need  not /be  used,    any  curve  of  the  correct  general  saape 
is  suitable  for  a  starting  guess,   and  the  better  the  guess,   the  more 
rapid  the  conversion  to  an  acceptable  solution.     Perhaps  the  best 
choice  would  be  fche  pattern  twWHk  a  similar  sh*p,   i±  such  is 
available. 

The  resulting  orthogonalized  pattern  will  now  be  used  as  a 
starting  point  of  the  iterative  process. 

(6)  Inertia  loading  uased  on  this  pattern  is  computed  and  plotted    -   -  a 
function  ot  X,  L  along  the  length  of  the  ship. 

(7)  This  is  integrated  to  obtain  shear,    to  which  is  added  the  routory 
inertia  term,   giving  the  derivative  of  bending  moment.     (In  the 
first  cycle  of  each  mode  shear  is  the  diiiereatial  ol  bending 
moment).     M'  is  plotted. 
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($)    M*  is  integrated*   and  the  integral  curve  is  divided,   at  each 
station,     y  flexural  rigidity  to  obtain  curvature,  which  is 
plotted. 

(9)  Curvature  is  integrated  to  obtain  the  bending  component  oi  slope 
to  within  a  constant. 

(10)  Shear  is  divided  by   shear  rigidity,  and  the  resulting  shear  slope 
is  added  to  the  Lending  *»aog  oucnt,   and  the  total  is  used  in  an 
auvi  roces*  to  determine  a  constant  whicu  locates  the 
axis  (the  constant  oi  the  integration  in  (9)),  which  is  added  also, 
to  give  total  slope  .     This  is  plotted. 

(H)    Slope  is  integrated  to  give  the  deflection  pattern  *.o  within  an 

additive  constant,  which  is  used  in  an  auxiliary  process  to  obtain 
the  axii-  shift.     The  resultant  defection  pattern  is  tabulated. 

(12)  I  his  pattern  is  now  submitted  to  the  same  orthogonality  checks 
used  on  the  starting  guess.     Ihe  corrected  curve  is  normalised 
to  give  a  pattern  whose  value  is  unity  at  the  center  -  or  center- 
most  -aiter-antinode.     The  unnormalixed  value  at  this  antioode 
is  retained  lor  use  in  the  xrequency  calculation.     Ihe  normalised 
pattern  becomes  the  sasi*.  oi  the  inertia  loading  Lor  the  next 
cycle,   or  is  at.  d  as  the  final  pattern,      ibis  decision  is  based 
on  the  amount  oi  change  oi  deflection  pattern  and  frequency  be- 
tween this  cycie  and  the  previous  one. 

(13)  Frequency  is  obtained     .   equation  (2),   where  (  y*)»tart  s    *»   *°d 
(y*)  final  is  the  unaormalized  antinode  deflection  retained  in  (12). 

(14)  In  the  flail  cycle  of  each  mode,    the  curvature  pattern  i« 
normalized  by  dividing  throughout  by  (y*)  final  and  retained  as 
the  final  curvature  pattern. 

(15)  The  shear  slope  distribution  is  subtracted  from  the  corrected 
total  slope  distribution  to  give  the  corrected  bending  slope,     ihis 


«i 


mm 

Si 


»q«ia  to 


l 


tdw  |« 

.-.  (- 

Iftta*". 

- 


■fl 


•i  ^  gaisa#4  £■ 


CI) 


,20 


is    normalized  b-y  division  ay  (y*)  final  and  retained  for 
use  in  the  rotatory  inertia  correction  of  the  subsequent 
cycle. 
(16)  The  pattern  obtained  in  (12)  is  now  used  to  repeat  the 
entire  cycle*  if  deemed  necessary. 

pc-ndix  A,    Part  1  contains  the  derivations  of  the  important  relations 
used.     Part  2  is  a  detailed  outline  oi  the  procedure,  containing  formulations 
of  all  operations  as  used  in  the  application  of  the  method  illustrated  in 
Appendices  5  and  C. 

ihe  method  was  executed  using  structural  data  ior  a  three-ei&idhs 
scale  model  of  a  sub       rin«,   and  three  lateral  mode  a  were  computed*   each 
having  ueen  carried  through  three  cycles  of  iterations,   this  being  sufficient 
to  show  convergence  oi  the  method  and  to  reduce  residual  error,  as 
indicated  roughly  by  the  last  increment  oi  change,   to  an  acceptable  level. 


>:, 


•idJ  fc«fliil«-r   wi* 


mn  m%»i  m  *" 


j?  ,.1    ij  ■   *      '•  «•■ 


Y"!  »*<?■- 


<    ia  ■.*>-!  .'■■  Ml i   :"<.'.  '    hfl    1JK'.  t  • 


.£-  *a« ■■• 


^•m*  ulu  I 


i    ■ 


. 


«*i  Ml 


aa 


in      *k.b\ji,i 


.  y    observation  oi  the  convergence  of  the  results  oi  iteration, 
patterns,  normalised  to  give  &  value  of  unity  at  the  center     or 
centermoet   after   antinode,  are  presented  in  Figure  XXXI11.     The 
corresponding  ^u- valuta  .attexos,   froa.  which  strain  patterns  can 
he  outlined,   were  normalised  in  such  a  manner  that  when  the  factor 
oetween  normalised  deflection  and  true  deflection  wnarta—  has  ceen 
obtained,   this  same  fas  tor  multiplied  by  the  normalised  curvature  ■tWf 
will  give  the  maxima  of  the  true  curvature  pattern,   a.i  di     ensions 
being  consistent.     These  patterns  are  presented  in  Figure  XXXIV. 

Table  I  gives  the  final  f requeue ies  and  periods  of  the  three  mode* 
vaiculated  for  the    node 

HATftfAc,  TBANSV  F*n.QU£Nwla.i»  JhUL>  .i-iJSlOOS 

WQ  &.  S  FOOT  SUBMS.it QM.D  St 


Natural  F  req. 

1 

h *UO 

riod  (  milliseconds) 

417 

1 

143.9 

Mi 

2.01 

71.5 

13V1  • 

3.  34 

43.  1 

I 
I 
y 

•    See  supplementary  discussion,    App«ndi>.  D 

Table  II  tists  the  offsets  of  the  normalized  deflection  and  curvature 
patterns  pictured  in  Figures  XXXIII  and  XXXIV. 

s  to  the  method  developed,   it  was  found  to  ce  convergent,   though 
it  appears  that  the  delayed  introduction  of  rotatory  inertia  into  the  pro 
NM  must  act  to  delay    convergence.     Thi*  is  compensated  for  by  the 
frequent  orthogonality  checks,    especially  their  application  to  th. 
starting  pattern. 

The  results  and  method  are  discussed  in  the  subsequent  section. 
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iMMtum^mssmaL 


x/h 


1 1 

*  1 


tt 


y> 


2 


MMHW 


it 


0,000 

0,025 
0*050 

0.075 
0.100 

0.125 

0.150 

0.175 
0.200 
0.225 

0.250 

0.275 
0.300 

0.325 
0.350 

0.375 
5.400 
0.425 
0.450 
0.475 
0.500 

0.525 
0.550 

0.575 
0.600 

0.625 

0.650 

0.675 
0.700 

0.7*5 
0.7J0 
0.775 
0.^00 

0.825 
0.850 

0.( 
0#< 

0.950 

0.975 
1.000 


-3*408 
3*0; 
2.6* 

1.994 
1.  ">#1 
1. 

0.1 
0*682 

0.395 
-0.125 

+0.109 
0.325 
0.511 
0.661 

o,7aa 

0.878 
0.944 

0.987 

1.000 
0.9*9 

0.965 
0.912 
O.835 
0.711 

0.56$ 

0.3§3 
+0,169 

-O.Ofc? 

0*350 
0.454 
0.976 
1.310 
1.662 

2. 3 

2.783 
3.207 
3.624 
4*066 
•*4*502 


0.' 
0.81 
0.706 
0.532 

0.336 

-0.133 

+0.095 

0.302 

0.51> 
0.69 3 

0.83a 
0.93s 

0.995 

0.987 

0.9&* 
0.780 

O.582 

+0.304 

-0.024 

0.411 

0.849 

1#227 

1.8' 

3.808 
-4.45C 


V3.925 
3.116 

0.729 
+0.102 
-0.417 
0,806 
1.074 
1.178 
1.150 
1.071 
0.862 

0.596 

-0.295 

+0.036 

0.330 

0.602 

0.793 
0.928 

0.996 
0.975 
0.875 
0.677 
0.417 
40.083 

-0.275 
O.655 
1.002 
1.227 
1.373 
1.3«9 
1.281 

0.951 

-0.487 

+0.155 

1.044 

2.297 

3.711 

5.154 

+6.590 


' 


. 


23 


25 

IV.     DISCUSSION    OF     .    -       jj  L  i  b 

A*     Mode  Pa  tic   ns: 

The  shape  of  the  mode  ^ttcrai  (Figure  AX XIII)  gtner  tt  the. 

feutho    *'    — ycllll— i.     In  ships  and  »u  mes,    uoth  the  OBftM  and 

fie*u-&l  rigidity  are  greater  at  the  midsection  than  at  the  ends.     This 
has  the  effect  of  altering  the  mode  shape  so  that  the  center  deflections 
a*e  smaller  relative  to  the  end  deflections  than  in  a  uniform  beam, 
and  the  nodai  points  are  closer  to  the  midship  section.     For  example, 
for  the  first  mode,    the  venter  deflection    naxi  .  urn  of  the  uniform  beam 
amounts  to  60 %  of  the  end  deflection,   while  in  the  esse  of  the  submarine 
studied,   the  anti  node  deflection  was  29.4'/«  of  the  bow  deflection  and 
ZZ.H%  of  the  stern  deflection.     The  nodes  for  the  uniform  oeam  otvur 
at  ZZ'k  from  the  ends,   while  for  the  submarine  investigated,   they 
occurred  at  Zb.Zb  from  the  oow  and  30.  7%  from  the  stem  respectively. 

Figure  I  clearly  demonstrates  the  proncuxed  increase  in  mass, 

sectional  area,    moment  of  inertia,   and  shear  factor  ft  in  the  vUinity  of 

the  hard  tank  and  conning  tower  amidships  in  the  case  of  the  submarine. 

One  can  readily  arrive  at  two  consequences  of  this  discontinuity  of 

structural  strength.     The  first  of  these  has  already  been  mentioned  - 

that  of  the  ve   y  high  end  deflections  relative  to  those  in  the  middle  tody. 

This  behavior  will  obviously  occur  for  al;  the  modes,  with  the  resulting 

trend  towards  patterns  with  many  small  loops  in  the  middle  with  large 

end  branches.     It  is  easy  to  see  how  any  loading,   such  as  an  explosion, 

M    urring  in  the  middle  body  of  a  submarine  can  incur  very  large  deflec 

tions  at  the  ends  with  attendant  small  deflection  near  the  load.     In  the 

case  of  the  third  mode,   the  deflection  at  the  stern  will  ue  six  Uxorcs  that 

at  the  center  I    The  reader  will  readily  appreciate  the  importance    ° 
this  characteristic  of  the  mode  patterns  with  respect  to  the  bodiiy 
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response  of  fuo^mrioei.   Severe  whipping  :v,ay  be  expected  to  occur 
whenever  the  contributions  oi  the  odd  end,  or  even  numbered  mode  end 
deflections  become  additive  for  any  general  vibratory  emotion  of  the 
submarine. 

The  second  consequent  of  the  irregular  distribution  of  structural 
characteristics  becomes  manifest  in  the  effect  of  rotatory  inettia    and 
shear  flexibility.     i  or  a  given  structure  utass  per  unit  length,   a  sub- 
marine will  he  maay  timet  stiffcr  than  a  solid  prismatic  oar,   since 
in  the  case  of  the  former,   the  crossection  is  essentially  a  circular  ring 
of  material  relatively  distant  fro  n  its  centroid.     This  fact,   coupled  with 
that  of  the  nature  of  the  longitudinal  distribution  of  fiexural  rigidity, 
explains  why  the  .magnitude  of  the  rotatory  inertia  correction  to  the 
frequency  of  the  first  mode  is  somewhat  greater  for  the  submarine  than 
for  the  uoifor  xi  beam  or  even  a  surface  ship.     Its  increased  moment  of 
inertia  for  a  given  mass  ft  is  compared  to  the  uniform  beam  results  in 
larger  values  of     im     .     £ven  for  equal  slopes,   the  rotatory  inertia 
correction  would  be  larger  for  the  suiu  .  arine.     However,   the   in- 

creased flexibility  of  the  end  sections  over  their  uniforu    beam  counter- 
parts causes  the  ttof  t 9$  increase  over  all  the  length  of  the  sumr.arine, 
the  greatest  increase  occurring  between  the  ends  and  the  nodes.     The 
steepness  of  all  the  deflection  patterns  near  the  ends  results  in  a  rather 
selective  effect  of  rotatory  inertia  upon  the  patterns.     In  general,   the 
effect  of  rotatory  inertia  upon  the  mode  patterns  studied  was  most  pro- 
nounced at  the  ends,   and  resulted  in  an  increase  in  the  relative  deflection* 
ison  of  experimental  and  theoretical  mode  patterns  for  a  similar 
submarine  where  rotatory  inertia  was  neglected  in  the  latter  calculation, 
shows  larger  experimental  end  deflections  than  the  theory  predicted. 
The  authors  feel  that  such  a  discrepancy  may  te  due  to  the  application 
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of  an  incomplete  theory,   particulaty  in  the  light  of  the  results  of  this 
thesis. 

Shear  flexibility  apparently  causes  little  change  in  the  shape  and 
magnitudes  of  the  normalised  mode  patterns.     The  shape  of  the  KAG 
curve  (Figure  I)  indicates  that  the  location  of  maximum  shear  near 
the  quarter  points  of  the  submarine  will  produce  the  largest  value  of 
shear  slope.     Actually  the  percent  contribution  to  y*  of  shear  slope 
was    nu.h  more  uniform  over  the  full  length  of  the  submarine  than  was 
expected.     This  accounts  for  the  small  effect  of  shear  deflection  on 
the  normalised  patterns.    As  was  expected*   the  value  of    (3       increases 
with  increasing  modes  numbers,   and  thus  the  effect  of  shear  deflection 
ttpon  frequency  becomes  more  important  with  increasing  mode  number. 
frequency  behavior  is  discussed  more  fully  in  a  subsequent  section. 
It  should    be  mentioned,  however,   that  with  the  third  mode,   the  percen- 
tage of  shear  component  of  slope  gave  indications  of  Lecoming  larger 
in  the  region  of  the  middle  oody,  where  the  greatest  effect  on    cJ*  will 
result. 

**•     Curvature  Patterns; 

The  curvature  patterns,   Figure  XaXIV,   clearly  indicate  possible 

regions  of  structural  weakness  with  respect  to  resistance  to  damage  from 

whipping.     In  particular,   it  appears  that  excessive  strain  can  be  expected 

at  scout    X/L    *    0.9  in  the  event  of  third  mode  whipping  of  significant 

magnitude*    The  magnitude  of  the  strain  wire,   of  course,   also  depends 
upon  the  radius   of  the  su&marinc  hull  at  this  location. 

*    Note  that  the  rudder,    propellers  and  shafting  add  mass  with  little 
attendant  increase  in  stiffness  in  this  part  of  the  sut marine. 
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Except  for  the  apparent  deficiency  just  noted,   the  structural  design 
of  the  suo marine  appears  to  ce  quite  well  balanced.     The  general  shape 
of  the  curvature  patterns  obtained  by  the  authors  agrees  quite  favorably 
with  experiments -ly  determined  curvature  patterns  of  similar  submarines. 

C.  r" requeue ics  : 

Frequency  oehavlor  observed  in  this  investigation  generally  confirmed 
expectations  Lased  on  theoretical  consideration  of  the  secondary  elastic 
effects.     Lxact  breakdown  of  these  effects  is  not  available  in  that  no  com- 
plete iterative  calculation  was  conducted  without  these  effects.     Although 
the  first  cycle  of  each    node  involved  no  rotatory  inertia  effect,   it  cannot 
be  considered  as  yielding  a  useful  figure  for  comparisons,    because  of  the 
predominating  effect  of  the  starting  estimate  curve.     However,   if  we  accept 
these  first  cycle  frequencies  for  a  qualitative  examination  (Table  111, 
Appendix  l),   we  see  the  expected  drop  in  frequency  when  rotatory  inertia 
is  introduced,   of  magnitude  greater  than  what  can  be    attributed  to  con- 
vergence alone.     It  further  appears  to  be  of  increasing  importance  with 
higher  modes. 

Dsta  showing  the  independent  effect  of  shear  flexibility  is  not  available, 
but  preliminary  investigations  summarized  in  reference      [?;]     showed  a 
similar  depression  in  frequency  when  this  effect  was  introduced,   although 
the  effect  was  separated  for  the  first  mode  only.     This  is  generally  con- 
sidered to  be  of  increasing  importance  for  higher  modes. 

The  cumulative  result  of  these  effects  is  noted  in  the  fact  that  exper- 

entally    determined  natural  frequencies  correspond  more  nearly  to  the 

1       ?       2 

scries  1,  2,  3,   than  to  the  3**.    5*,    7   ,   scries  obtained  theoretically  from 

a  uniform  beam  without  secondary  effects.     fi?J       This  same  tendency  is 
o  .served  in  the  frequencies  obtained  in  this  investigation,   which    were  seen 
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in  Table  1  to  correspond  to  the  series  1,   2.  ol,    3.34,   or,   for  contrast, 

2  2  2 

expressed  in  the  second  form,,     3   ,   4.25   .    5.2   . 

However,   comparison  with  results  of  earlier  computations  for  the 
same  model,   but  without  the  rotatory  inertia  effect,   cast  some  douct  on 
the  accuracy  of  the  third    node  frequency  of  1 3  > i  Ipsa,   indicating  it  is 
perhaps  too  high.     This  suspicion  is  reinforced  by  the  3.  34  value  in  the 
above  series.     The  error  (pro'jaoly  amounting  to  less  than  4%)  is  con- 
sidered to  be  due  to  the  fact  that  the  inter/node  orthogonality  check  as 
formulated  is  inexact  for  other  than  the  trivial  modes,   and  the  error 
becomes  of  some  importance  for  the  third  and  higher  modes.     The  more 
correct  formulation  is  developed  in  appendix  D,   and  the  third  mode 
frequency  corrected  by  its  partial  application. 

D.    Sources  of  *.;  ror  ; 

Several  sources  of  error  arise  from  fields  the  detailed  consideration 
of  which  is  beyond  the  scope  of  this  thesis,    out  are  noted  for  completeness. 

The  first  of  these  stem  from  the  information  included  as  Structural 
Data,  found  in  appendix  i>.     A  submarine,    though  more  regular  in  design 
than  a  surface  ship,   is  a  complex  structure,   and  the  participation  of  its 
components  in  carrying  the  loading  are  more  difficult  to  analyse  than 
simple  beam  theory  suggests.     Therefore  agreement  is  hard  to  find  as 
to  which  memuers  aiouid  be  included  in  calculating  moment  of  inertia, 
v.  ross  sectional  area  resisting  shear,   and  the  shear  distribution  eoeffi- 
went  K. 

Far  more  controversial  is  the  matter  of  entrained  water,    here  it  is 
especially  difficult  to  verify  a  given  theory  because  vibration  theory  is 
not  sufficiently  refined  to  permit  separating  out  other  effects  which  inter 
act  with  the  manifestations  of  entrained  water. 

.  oth  of  the  above  questions  have  been  sidestepped  ty  utilizing  data 
already  calculated  from  plans,   based  on  the    assumptions  outlined  in 


.i.«l«o»  «1  .t.  .M  .£   M  .t   .1  •»"»•  •*»  •»  *••*••'  ql 

*  aili 

-    .1    .aI  >i-   mum    li 

p  .1  (**  a«<i)  ••»!  of  s»ll0t  i<3# 


•Tt 


b*UA3»c  »d4  **>J*il  me**  »««*  v  «#••  l« 


»*     l^*r 


«*!•»»  si 

4  »  fcxwi  rJ 


I 

i©  in« 


!30 


Appendix  E.     it  is  not  to  be  interred  that  these  matters  are  insignificant; 
they  are  vital  to  any  prediction  of  vibration,   and  are  the  subject  of 
considerable  study. 

It    has  :>een  assumed  that  deflections  are  small  and  that  the  pro- 
portional It wlH  have  not  ueen  exceeded.     Standard  values  of  Young's 
Modulus  aud  Poisson's  £atio  have  been  assumed.     Any  departure  fro 
these  assu  options    in  a  ship  tested  experimentally  can  be  expected  to 
.ause  discrepancies. 

The  second  group  of  errors  can  be  traced  to  the  procedure  and 
to  its  execution. 

First,   failure  to  carry  the  iteration  through  sufficient  cycles  leads 
to  an  error,     Ihe  residual  error  after  a  given  nun  u>ex  of  cycles  will  be 
large  if  (a)    the  starting  assumption  of  deflection  was  grossly  in  error, 
or  (.->)  if  orthogonality  checks  were  not  properly  and  regularly  applied. 
An  error  was  introduced  also  by  the  application  of  the  inexact  orthogono- 
nality  formulation. 

Next,   the  rotatory  inertia  term,    weing  cased  necessarily  on  data 
from  a  previous  cycle,   is  slightly  in  error. 

Further,   the  failure  to  make  this  correction  in  the  first  cycle  re- 
sults in  a  discrepancy  which  will  take  several  cycles  to  eliminate  entirely. 

Shear  lag,   and  probably  other  lesser  effects,   have  not  been 
considered. 

Xntegraph    operation  involves    some  small  error.     This  is  periodi- 
cally filtered  out  by  boundary  conditions  and  other  checks,    but  some 
residue  necessarily  remains.     In  particular,   a  small  amount  of  lost 
motion  was  noted  in  the  inte graph  used. 

Any  error  in  developing  a  given  mode  pattern  will  result  in  a  sir.ali 
error  in  each  subsequent  pattern,    through  the  interniode  orthogonality 
.hecks,   whKh  necessarily  employ  the  ordinate*  of  previously  determined 
nodes. 
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if  inaily,   the  usual  remarks  concerning  human  error  in  commuting 
values,   and  in  plotting,   integrating,   and  reading  curves.     It  is  noted, 
however,   that  the    ..ethod  by  nature  tends  to  reject  such  error; or  to 
make  them  immediately     bvious  as  the  technique  becomes  fa  miliar. 

£.  Applicability: 

The  method  of  this  paper  is  not  envisaged  as  replacing  the  rather 
advanced  techniques  of  digital  computers  and  electrical  analogs.     It 
may  well  coropleoMMSt  tnem  to  »an€  degree,   but  is  intended  primarily 
to  serve  the  needs  ol  |kt  smaller  research  facility  which,   operating  oa  a 
limited  budget,   is  faced  with  the  necessity  of  performing  such  determinations 
with  simpler  tools,     late  graphs  of  the  type  used  are  found  in  a  nun   er  of 
educational  institutions,   design  offices,   and  shipyards,    but  the  authors 
have  reason  to  relieve  that  more  than  a  few  of  them  lie  neglected  beneath 
at  least  a  decade's  accumulation  of  dust.     The  only  other  tool  needed  is  an 
automatic  or  semi-automatic  commercial  desk  calculator.     It  is  considered 
that  tne  work  can  bi  set  up  on  appropriate  forn  s  in  such  a  manner  that 
a  mathematics  aid  or  junior  research  assistant  can  carry  it  out  without 
special  training  and  onder  only  a  modicum  of  supervision. 
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V.       CONCLUSIONS 

Convergence    of  the  iterative  method,  with  secondary  efiects, 
as  developed  in  this  thesis  appears  satisfactory,   and  any  tendency 
towards  instability  is  held  in  check  by  periodic  application  of  orthogona- 
lity conditions.     The  method  inherently  detects  and  alters  out  errors, 
so  that  a  mistake  in  an  early  cycle  does  not  invalidate  subsequent  cycles 
of  the  same  mode,   «ut  may  delay  convergence  slightly. 

While  the  procedure  appears  complex,  it  can  ue  programmed 
on  standard  forms  for  application  oy  research  assistants  without  special 
knowledge  of  vibrations  or  advanced  mathematics.     Equipment  required 
is  standard  to  many  installations  and  inexpensive  relative  to  such  special 
equipment  as  digital  computers. 

Results  of  this  method  appear  to  be  acceptable;    the  third  mode 
frequency  is  considered  accurate  to  within  Z.6%  after  three  iterative 
cycles,  while  earlier  modes,  having  attained  better  convergence  in  the 
same  number  of  modes,  are  considerably  more  accurate. 

The  secondary  effects  are  considered  to  have  oeen  adequately 
accounted  lor,  with  little  added  complication  to  the  process,     la  particu- 
lar,  rotatory  inertia  appears  to  have  had  a  somewhat  greater  ef.ect  than 
has  been  evidenced  in  studies  of  uniform  beams;    end  deflections  are  seen 
to  be  increased  relative  to  antinodal  excursions  when  the  effect  is  con- 
sidered.    Secondary  eH'ects  thus  are  included  not  only  in  frequency 
determination,    but  also  in  deflection  and  curvature  patterns. 

clear  indication  of  highly  loaded  sections  is  presented  by  the 
normalized  curvature  distributions  which  result  directly  irom  the  pro* 
cess.     Structural  inadequacies  are  thus  pointed  out. 

It  thus  appears  that  the  method  developed  is  well  suited  to  the 
purposes  of  small  facilities  conducting  vibration  research  on  limited    ud- 
gets,  where  accurate  theoretical  determination  of  information  on  natural 
modes  is  needed. 
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VI      B  L  COMMEND  A  1  IONS 

It  is  recommended  that  the  method  for  modal  calculations  which 
is  described  in  this  thesis  be  considered  hf  small  research  activities 
for  similar  work.     As  a  general  rule,   such  activities  operate  under  a 
rather  limited  budget  with  correspondingly  small  allocations  for  co;n- 
puting  equipment  of  the  analog  or  digital  types.     On  the  other  hand,   most 
of  these  activities  will  certainly  own  more  than  one  desk  calculator  and 
will  have  the  funds  available  for  the   purchase  of  an  integraph  in  the  event 
that  one  is  not  owned.     At  least  for  the  modes  of  interest,   the  method  used 
by  the  authors  appears  more  than  satisfactory  from  the  standpoint  of 
moderate  equipment  costs,    simplicity,  and  availability.     Digital  compu- 
ters are  in  ever  increasing  demand,   and  the  delay  involved  in  obtaining 
results  from  a  central  calculating  contractor  may  be  unacceptable  for 
the  smaU  research  organisation.     The  method  is  clearly  not  competitive 
with  that  of  the  digital  computer  when  the  latter  is  readily  available  and 
a  large  number  of  mode  patterns  is  desired. 

Experimental  corroboration  of  the  authors'  results  would  >-<e 
highly  desirable     before  the  method  is  accepted  for  general  application. 
For  any  given  submarine  design  this  can  ~»e  accomplished  by  mechanical 
excitation  of  a  model  of  the  proposed  vessel,   provided  it  is  of  sufficiently 
large  scale  to  assure  structural  similarity. 

It  is  further  recommended  that  the  effect  of  rotatory  inertia  as 
well  as  shear  deflection  e>e  included  in  the  natural  mode  calculations 
regardless  of  the  ■  ethod  used. 

In  any    future  amplication  of  this  method,   it  is  recommended  that 
the  more  exact  orthogonality  conditions  which  arise  from  the  inclusion 
of  rotatory  inertia  effects  be  applied.     While  it  is  more  convenient  to 
negiect  them,   such  practice  may  result  in  the  loss  of  refinements,   espec- 
ially in  the  higher  modes,  which  are  of  the  same  order  of  magnitude  as 
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the  rotatory  inertia  corrections  themselves. 

The  effect  of  sheer  lag  ia  ship-type  structures  upon  their  vibration 
behavior  should  be  investigated. 
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AFPLNDIX         A 
PABT    1 

Theoretical  basis  : 

The  submarine  it  initially  considered  to  oehavc  as  a  long,   thin 
beam  with  vanishing  shear  force  and  tending  moment.     Under  this  acsump* 
tion,   then,   tne  classical  theory  of  the  lateral  vicratioa  of  a  free-free  non- 
uniform bar  applies,     C    3  T  ; J  •     We  are  concerned  here  with  the  normal 
modes  of  vi   ration.     When  a  bar  vi urate*  in  one  of  its  normal  modes,   the 
deflection  at  any  iocation    varies  harmonically  with  time.     If  damping  is 
neglected,   the  lateral  displacement  o£  the  oar  can  oe  represented  &y  : 

y      s    X   (    A     co*    *Jt    +        B     sin  Alt    )  (6) 

where  ^  is  a  function  of  the  coordinate  x,   and  determines  the  tfiape  or 
pattern  of  the  normal  or  natural  .node  of  vi;  ration  which  we  are  considering. 
)£     is  known  in  the  literature  as  a  "normal"  or  "natural"  function  . 

It  is  obvious  that  any  particular  normal  mode   ai  vibration  is  charac- 
terised solely  •  y  Jfj    and  the  frequency  of  that  mode,      «J^     ,     iry  super- 
imposing all  possible  normal  modes  of  vibration,   the  general  expression 
for  the  irec  lateral  vi oration  becomes: 

y     *      2        Xj    U|««t    *>i    '  *4sin*tt) 

/»* 

in  order  to  calculate  a  normal  mode  pattern,   one  must  resort  to 
beam  flexure  theory.     The  ample  (Juluier)  theory  of  sending  may  wt  assumed 
to  hold  if  the  vibration  occurs  in  one  of  the  principal  planes  of  figure     °* 
the  oar,   and  if  the  crossectional    dimensions  of  the  bar  are  small  compared 
to  its  length.     In  the  actual  case  of  a  ship  or  submarine  the  first  of  these 
restrictions  offers  no  problem,   but  the  second  can  M  unrealistic,   particularly 
for  a  modern  submarine  with  its  "cigar"form  and  essentially  circular  sections. 
Strictly  speaking  then,   the  i-.uler  wending  theory  is  not  applicable  to  a 
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submarine  or  even  a  beam  whose  cross  sectional  dimensions  are 
appreciable.     For  the  present*   however,  we  shall  use  the  kuler  theory 
as  a  beginning  and  modify  it  subsequently  by  a  correction  for  shearing 
force  and  rotatory  inertia  effects. 

We  refer  now  to  Figure  A-l,   which  is  a  representation  of  a  non- 
uniform beam  and  an    element  thereof.     The  symbols  are  in  accordance 
with  the  Nomenclature,   except  that    w    here  is  taken  as  some  loading 
in  lb/ ft.     With  the  axis  oriented  as  shown,   the  well  known  Eulcr 


differential  equation  for  the  deflection  curve  is: 

2. 


±1 


dx* 

FIGUR&  A-d 
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We  differentiate  (8)  once  to  obtain  the  shearing  force: 


dx 


(""&) 


»    -Q. 


An  additional  differentiation  yields  the  loading. 


dx 
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Equation  (10)  ie  the  differential    equation  of  a  oar  subjected  to 
a  continuous  load  of  Intensity    W   ,    In  the  case  of  free  lateral  vibration, 
one  can  apply  D'/Uembert's  principle  and  replace   W  by  the     inertia  load- 
ing 

m     fry  (11) 

O    t* 

Substituting  (li)  for  W  in  equation  (1-,),   the  general  equation  for 
the  later*!  vibration  of  the  jar  becomes 


Finally,    since,   for  simple  harmonic    notion.iy     =    -    yo)v  (13) 

we  may  write  equation  (12)  as  follows 


>  -£v).  M*  a* 


We  now  return  to  the  reality  that  the  dimensions  of  the  coss  sec- 
tion are  not  small  with  respect  to  the  length  of  the    *ar.     Investigation  of 
the  effect  which  the  cross  sectional  dimensions  can  have  on  the    frequency 
are  incompassed  in  the  classical  literature.     [63L^3C^5  J  .        Of  these 
perhaps  the  most  widely  accepted  is  the  treatment  of  the  effects  of  rotatory 
inertia  and  shearing  force  by  Timoshenko  in  reference    £.*•>  ~}L  55  ]  »    There 
are  douttiess  other  effects  introduced  by  the  cross  sectional  dimensions 
and  configuration  but  for  the  purposes  of  this  thesis  the  Timoehenko 
mechanism  was  adopted  sMtmM  it  accounts  for  tne  two  major  effects 
mentioued  above. 

The  effects  of  both  rotatory  inertia  and  shear  deflection  become 
of  considerable  importance  in  the    higher  modes  of  vi. ration  where  the  oar 
is  subdivided  into  shorter  segments  uy  an  ever  increasing  numotr  oi     -odes. 
[9   3C691 
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a.     Rotatory  Inertia  : 

Refer  again  to  figure  A-l.     One  can  easily  see  how  the  element 
dx  during  vibration  not  only  performs  a  translatory  motion,    but  can  rotate 
as  well.     The  variable  angle  of  rotation  is  equal  to  the  slope  of    the 
deflection  curve,     dy/£x,   &.nd  the  corresponding  angular  velocity  and 
angular  acceleration       will  be  given  by 

^2y      ,     and  ^3y 


The  moment  of  the  inertia  forces  of  the  element  aoout  the  axis 
through  its  centroid  perpendicular  to  the  x-y  plane  will  evidently   >e: 

/  JI     dx  \/     m       +        a  (ny  A)  dx    \f         tfy    \dx  (15) 

(The  moment  is  taken  positive  when  it  is  clockwise),   where  m/A  is  the  mass 
per  unit  volume  of  the    non-uniform  bar. 

We  now  expand  (15),    and  obtain: 


dMr 


1  m     (pe)  **  '  (dx)    l*T  2A  T    cy 

.     >I_    .       ^(TW/A)  1     \#±  (15  a) 


If  we  neglect  second  jrder  and  higher  infinite  si  i^als, 

Remembering     again  that  for  simple  harmonic  motion         — =      -  cS^Y*   we 

o     t 

finally  obtain: 

dM      -    -    Im        h    (-fc^y)**      =       Vm£    <*y  dx  (16) 

"a- "air  *      <^ 

This  moineat  should  be  accounted  for  in  considering  the  variation 
of  bending  moment  along   the  axis  of  the  bar.     Thus,    instead  of  equation  (4) 
we  will  have: 

d  M         ■       Q       +       I  mCJ*  hx_  (17) 
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Substituting  this  in  the  equation  to*  the  deflection  curv«, 


and  using  expression  (U)  and  (12),  we  finally  ©;  tain. 


*r&)    *»  "^  **■*  *  « 


b.     i.ffect  3  CS£ 

•U  irate  differential     equation  i*  obtained  II  not 

only  rotatory  iv  hear  dcflettioa  in  toosidered.      The  all 

of  the  del  -  curves  depends  eat  W  on  the  ms 

oi  the  tots  I  ut  M  Ihe  a^;  .  *«  due  to  shear  as  weii.       .et 

^     denote  trie  slo  >e  oi  the  vtf  when  shearing  fo±ce  is  neglected  and    ^ 
the  angle  of  shear  at  the  neutrAl  axis  of  the  »ar^e  rross  rection.     evi- 
dently th<  ihz  #um  of   'y    and^    ,   i.e., 

(20) 

iu»  the  eJ  i.ng  mo   seai 

and  rttoarfam  the  fallowing  expressions, 

(2i)    I    (22) 

is  a  nu  :*erit.ai  I  Is  a  fun.  of  the 

IM  se«.tional  area,   usu&ily  defined  a*  follows: 

*    average  shear  streaa  in  ar 

•ham*  stress  in  section  (23) 

ihis  is  the  static  shear  stress  distn   utioc  factor  us«*d   ,     Ti-r.oshenko  and 
discussed  in  atandard  texts. 
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The    complete  mechanism  of  the  bar  including  both    rotatory 
inertia  and  shear  deflection  effects  may  now  be  formulated,     hefcrring 
to  Figure  A-l  the  differential  equation  of  rotation  of  the  element  mam  n 

will  be: 

2 

*  dx       +     Q     dx     «     -i2^~  .   -4^-dx   •  (24) 


stituting  equations  (21)  and  (22)  we  ootain: 


(h-0 


2 


AG  I-m    b*  *    I 


(25) 


Neglecting  ell  but  first  order  infinitesimals,   the  differential  equation 
for  the  translatory  motion  of  the  same  element  in  the  horizontal  plane 
will  be: 

z 


or 


J&JL-     .         £_!«*•#£    -^L    0 


»-Ttr-  -     §r-Jr     V7i      '  ■  °-  <*7> 


Taken  together,   equations  (25)  and  (27)  compete  the  simultaneous 
differential  equations  which  describe  the  more  complete  mechanism  of 
lateral  vibration  for  non-uniform  osrs. 
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'i  imoshenko's  analysis  of  rotatory  inertia  and  shear  deflec 
tion  effects  for  a  uniioi  ...  bftJ  may  Lt  found  la  reference    L^J    • 
He  obtains  two  si ^  ultaneous  equations  similar  to  (25)  and  (27). 
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He  then  proceeds  to  eliminate  from  equations  (2c)  and  (29)    ^ 
and  obtains  there L>y  the  weii  known  fourtn  order  coupled  equation: 

m-Jul.   4       J!X  .  /Le^L  ♦  JS     \  ft  .     ♦ 


t .         ** 


V 

* 


K    AVG  J     '  (&) 


*    Tixnoshenko's  notation  ha©  teen  altered  to  conform  with  the  authors'. 
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liquations  (Zj)  and  (27)  appear  not  to  ue  susceptive  to  such  an 
operation*   since  I,    A,   m,  and  K  are  all  functions  of  x  in  this  case. 
Considerable  effort  was  spent  by  the  authors  in  an  attempt  to  eliminate 
y  from  equations  (25)  and  (27),   without  success.     If  the  operation  is 
possible,  and  this  the  authors  doubt,   the  resulting  equation  will  at  best 
be  too  co  and  un  wieldly      for  practical  application. 

The  difficulties  involved  in  the  solution  of  equation    (30)  are 
well  known.     It  is  quite  obvious  that  the  irregular  variation  of  I,    «,   K, 
and  A  with  X  in  the  case  of  a  ship  dictates  that  equations  (25)  and  (27) 
must  oe  solved  in  an  approximate  .uanner  usin^  the  graphical,   lumped 
parameter,   or  fcayleigb  rocedure.     1  he  task  is  so  formidable 

that  the  authors  could  find  no  exact  solution  in  tac  literature  encompas- 
sing both  shear  and  rotatory  inertia  effects  applied  to  an  irregular  non- 
uniform far. 

To  recapitulate,   the  oasic  equations  used  by  the  authors  are, 
in  simplified  uir.n: 

u*  y  (3i) 

Q    *  0     -     KACI    (y*       -    r  )  (32) 

M'  .    Q    +    _!"»'*  (33) 

r   .  -£-  (34) 

r  •  r  +  /*  us) 

The  boundary  .traditions  applicable  to  the  quadruple  integration  mm: 

Q    *    0  for  *    0,        x     ~  (36) 

U  «    0  for  *    0.        X     «  (37) 
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Equation*  (31)  -  (3?)  describe  an  eigenvalue  proolern,  with  toeing 
the  eigenvalue.     Solutions  will  exist  only  for  a  discrete  set  of   fc)v  .     For 
each  such  value  of  &>*~  a  unique  mode  pattern  will  correspond.     An  excep- 
tion to  this  is  the  case  when    6iv  »    0,  where  two  possiole  modes  exist: 

y    «    1    (rigid  body  translation),   and    y    »    1  -    (rigid  -jody  rotation). 

x 

The  two  modes  for  this  tunique  case  have  L*eeo  orthogooaiited  for 
convenience  . 

Useful  orthogonality  relations  exist  between  the  normal  mode 
patterns,     in  particular,   for  each 

/,-.  dx      ■      |  (36) 

/c 
x  Hi  y,  dx    »  (39) 

?- 

and  approximately,  for  the  lower  modes, 

L. 

fax  y.  y#  dx      «    0,  i  &  i  (40  )• 

Murray    {,111     Y  application  of  integration  vy  parts  to  equations 
(38)  and  (39)  discovered  another  interesting  orthogonality  relation  .    which 
has  apparently  not  oeen  noted  previously  in  the  literature. 


fl\     dx  f*  (     $       -       A)         df 


Hi) 


Equation  (41)  is  used  oy  Murray  and  the  authors  to  evaluate  the  constant 
of  integration  of  jVdx,  and  thhs  exactly  determine  the  location  of  the 
true  slope  axis  (This  axis  has  previously  oeen  located  arbitrarily  at  the 
mid-length  point  of  the  slope  curve). 


•    The  reader  in  referred  to  appendix  D  for  the  derivation  of  the  exact  inter- 
mode  orthogonality  conditions  and  for  a  discussion    as  to  the  applicability 
of  the  approximate  condition,    Equation  4.. 
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If  we  take    y    a    1    +      C     and  substitute  in  equation  (41),   it  is 

iay  be  < 


easily  seen  that  C„  may  be  evaluated  as 

3         rL  rX 


3 


*  An  (  ^    -      x)      df    dx    '  (42) 


tt 


In  a  like  manner,    use  is  made  of  equation  (38)  to  determine  C    ,   where 

y  -   1    *    c4. 

It  must  be  emphasized  that  orthogonality  among  all  the  modes  must 
be  rigorously  maintained.     l_2  J  |J71  J     A  starting  pattern  for  the  iterative 
cycle  is  not  strictly  admissable  unless  it  has  been  explicitly  adjusted  to 
conform  with  the  orthogonality  conditions  described  in  equations  (38), 
(39),   and  (40).     The  resultant  pattern  is  not  valid  until  the  same  adjust- 
ment procedure  has  oeen  rigorously  executed. 

Orthogonality  with  respect  to  rigid  body  translation  is  insured  by 
use  of  a  zero  shift  of  the  y-axis:      y.    =    y.    +    A.     It  is  clear  from 
equation  (38)  that:  L 

A   /  mdx        ■      -     /         y    m    dx  •  (^3) 


i   I  mdx        ■      -    /        "y 

I  1     ' 


Orthogonality  with  respect  to  rigid  body  rotation  is  insured  by 
use  of  a  rotation  of  the  y-axis:         y.      ■  +      £      (x-x). 


From  equations  (38)  and  (39)  we  may  write 

y.   (x      -     *x)    dx      ■      0,  (44) 


t 


m 


and  substitution  of  the  above  expression  in  equation  (44)  yields: 


T 

x)JH  dx      s  — |my.     (x      -      x)      dx  .  (45) 


-Ly. 


to  so 

It  is  further  required  that  ooth  the  starting  and  resultant  deflection 

patterns  be  orthogonal  with  respect  to  each  flexural  mode  prior  to  the 
one  of  immediate  interest,    we  formulate  the  expression 
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y.. 


'.  «   y.  ♦    A    ♦    B  <x-x)   +    C^  ♦    C2y2  +     6^,7,..,  (**>) 

Constant*  C,»    C     etc  are  readily  determined  ay  substituting  y=    •    y    +  Caj.i 
etc.   into  equation  (4d).     For  example,   if  we  are  interested  in  adjusting 
the  second  mode  so  that  it  is  truly  normal  to  the  first  mode  we  set 

y   *    y   ♦   C    y 
whence  from  (4 

CJm   y/dx     *     •  Jm  ft    ff   dx  (47) 


in  general 


\  [m  jr?*    ■  /       y.  v.   dx  (4ft) 

J*       }  J.      »   ; 


Not  until  the  ove rations  indicated  Uy  equation  (4£)  are  carried 
out  should  a  cycle  be  begun  nor  a  resultant  deflection  curve  tpted. 

Since  the  frequency  depends  upou  the  value  oi  the  resultant  curve  at  the 
normalizing  ?oint,   y.  *,   the  frequency  showed  be  calculated         ,    .*ftcr 
ortnogonau^ation  is  completed. 
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Part  II 

Detailed  Outline  of  k  rocedure  : 

This  section  supplements  the  general  outline  of  the  procedure 
of  this  thesis  presented  in  Section  II,   and  sets  down  the  mechanics  of 
the  operations  implementing  the  relationships  derived  in  the  theoretical 
section.    Part  1  of  this  Appendix. 

Some  remarks  helpful  in  interpreting  the  analytical  expressions 
and  the  curves  of  ..ppendix  B  arc  in  order,     ftoman  numerals     indicate 
figures  in  Appendix  ft, 

(1)  Ail  curves  plotted  from  data  given  or  calculated  are  plotted 
at  intervals  of  0.  02$  U  and  faired  through  all  spots.     This  is 
considered  to  give  the  least  number  of  stations  which  will 
give  a  sufficiently  accurate  determination  of  true  curve  shape. 

(2)  Since  the  length  coordinate  has  been  ndn-dimensiouaiized 
(»x,  JL),   each  integration  performed  with  respect  to  it  results 
in  a  factor  ot  L  in  the  denominator  of  the  integral  expression. 
These  factors  are  removed  at  two  points  in  the  cycle:    when 
determining  curvature,  and  when  normalizing  deflection.  How- 
ever,  their  presence  within  the  cycle  results  in  appropriate 
modification  of  correction  terms  calculated  to  add  to  functions 
in  the  cycle. 

(3)  Since  the  virtual  mass  distribution  of  the  starting  data  was 
in  pound  units,   the    m    of  the  theoretical  relationships  has  oecn 
replaced  oy    ma    ,        the    a     (acceleration  of  gravity)  has  been 
carried  throughout  as  a  constant  and  is  removed  at  the  normaliz- 
ing stage  and  in  the  frequency  determination.     The  a  constant 
necessarily  appears  also  in  the  various  auxiliary  functions 
involving      m    . 
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(4)  The  horizontal  backstepping  of  the  various  integral  curves 
in  the  figures  of  Appendix  B  is  inherent  in  the  operation  of  the 
mechanical  Integraph  used.     (See  Plate  I)    For  convenience,   the 
backstep  has  keen  set  at  integral  numbers  of  inches,   either  two, 
three,  or  four  inches,  and  is  readily  apparent  from  the  y-axis 
labelled  for  each  integral  curve.     Coordinates  along  the  horiz- 
ontal axis  have  not  been  reuumoered  for  each  axis  shift,    but 
should  be  transposed  mentally  in  examining  integral  curves. 
Vertical  shifts  of  integral  curves  are  made  where  presentation 
is  facilitated  thereby. 

(5)  Also  inherent  in  the  integraph  is  a  scale  reduction  of  the 
integral  curve  with  respect  to  its  argument.     Scale  factors  used 
in  this  investigation  ranged  from  Uaee    to  7.5.    while  not  indicated 
directly  on  the  curves,   they  oi  course  determine  the  scales  marked 
on  the  integral  axes.     The  continuous  maintenance  of  correct 
scales  is  vital  to  accurate  determination  of  frequency  from  the 
deflection  normalizing  factor. 

(6)  Figures  illustrating  the  method  have  been  so  arranged  that 
the  four  principal  integrations  comprising  a  single  cycle  appear 
in  one  figure  (e.g.  Figure  VI)  while  the  auxiliary  curves  per- 
taining to  the  same  cycle  appear  on  one  or  two  separate  sheets 
following  (e.g.  Figure  Vil).     At  the  end  of  each  mode  is  a  figure 
(  XII,   XX,   XXXII)  illustrating  the  convergence  of  the  normalised 
deflection  patterns  from  tiu.  orthogonali«ed  uniform  beam  pattern 
ultimately  deemed  acceptaole  as  the  normal  mode  pattern. 

Sequence  of  Operations: 

When  the  taoles  and  smooth  curves  (1)  of  structural  information 
have  been  prepared,  it  is  necessary  to  obtain  the  coordinate  of  the  longi- 
tudinal ceatroid  of  the    appropriate  virtual  mass.     This  may  t>e  done  by 
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any  standard  technique,    vut  figure  II  illustrates  a  simple  and  unique 
method.     The  tangent  to  the  end  of  the  second  integral  of  the  mass  distri- 
bution is  extended  to  the  usee  line*   and  the  resulting  intersection,    rela- 
tive to  the  origin  of  the  second  integral  curve,  is  the  abscissa  of  the 
longitudinal  centroid  of  the  mass.     The  proof  of  this  is  given  in       f  7]  . 

Sufficient  information  is  now  on  hand  to  prepare  the  tables  of 
working  functions,  including 

mg     ,  mq  l/LTA,      10   L  /gfc   KAG/10   ,     (£-    £)/ L. 

rK  2 

g  (  £     x)/L.   and        gjm  (  £   -x)dJ/L 
As  they  become  known,   there  will  be  added  to  the  working  functions 
table  the  following  tabulated  functions : 

Yi     Ylrn     *     y2'     Y2  m        '   *tC* 
This  table  is  prepared  for  intervals  of  0.025  L. 

Certain  definite  integrals  and  other  constants  useful  throughout 
the  application  of  the  method  are  also  calculated  and  listed  at  this  point: 

o  dx/L  (49) 

rL  2  3 

N      «         J  (x    -    5)    mtf  dx/JL  (50) 


fjKi-      *)    mg  dy    dx/i.4  (51) 

Subsequently,  when  data  is  available,   the  following  met    added: 


w,t    I y    m 
1      J4    i 


o-dx/L  (52) 


2 

m<f  dx/i-,   etc.  (53) 


it  is  now  possible  to  commence  the  process.     While  any 
approximation  to  the  mode  shape  can  be  used  as  the  starting  curve,   it 
is  found  that  the  most  accurate  approximation  to  the  final  mode  shape 
available  will  lead  to  accelerated  convergence,   hence  fewer  cycles  of 
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iteration  for  the  de aired  accuracy. 

The  proposal  here  is  to  take  a  known  pattern,   as  from  a  similar 
ship  or  a  uniform  beam,   and  force  it  to  comply  with  conditions  of  ortno- 
gonality  applicable  to  the  ultimate  normal  pattern.     Thus  in  the  investiga- 
tion conducted  for  this  thesis,   the  normalized  ordinatss  of  the  patterns 
found  fcf  ...ord  iayleigh  are  taken  arbitrarily  as  a  starting  point. 

The  following  description  traces  the  development  of  a  tmrd,   or 
four-noded  mod*,   to  demonstrate  the  typical  case  of  a  higher  mode. 

From  the  working  function  tacle  and  the  values  of  Y    ,   tabula- 
tions are  obtained  and  plotted  oi  '"??••    mg(x  -  x)  Y   /  L,    mgy,    Y    ,   and 
ro«y?Y    .     i  he  four  resulting  curves  are  integrated  (aaII),   and  the  values 
obtained  are  divided  by  the  operational  constants  indicated  to  obtain 

A      ,       .      A    "*    »»  *NU  (54) 
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(56) 


,    "X  mV*  Y»      d*   ~ (57) 

ii  y     ' 


Now  the  values  of  (x    -    x),  L  are  multiplied  bf      B      the  values 
of  y    are  multiplied  and  the  values  of  y     are  multiplied  oy  C    .     Then 


obtain 

y      •     A      +      &      (x    -    x)  /  L     ♦      C.y       +      CtiTI 

All  values  of  y  are  divided  oy  y*,   its  value  at  its  normalizing  point  (the 
center -or  centermost-after  antinode  as  obtained  ty  i  plot  of  that  region) 
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to  obtain  a  tabulation  of  y  .   the  orthogonali2.ed,   normalised  starting 
curve  for  the  first  eycle    of  iteration.     (But  see  discussion  of  orthogo- 
nality in  Appendix  D). 

The  cycle  itself  (XX ill)  commences  with  the  curve  of  inertia 
loading    «      ■    ymo     .  which  is  calculated,   plotted,   and  integrated,   to 
give  a  function  of  shear, 

gQ/i-      »  -  y   ag-  (59) 

where  the  constant  oi  integration  is  &ero,   since  the  integral  curve 
starts  at  the  origin,   and  the  boundary  condition  of  *ero  end  shear  dictates 
that  the  true  shear  curve  should  also  pass  through  the  origin.     However, 
the  right  end  oi  the  integral  curve  may  differ  from  aero  by  a  small 
amount  (usually  less  than  5%  of  the  maximum  ordinate  of  the  curve). 
When  this  occurs,   a  linear  correction  is  applied  to  the  curve  to  again 
conform  to  the  above  boundary  condition.     This  is  done  by  drawing  a  new 
axis  joining  the  ends  of  the  shear  curve,   measuring  its  ordinate*  from 
this  curve  out  normal  to  the  old  axis,   then  replotting  on  the  old  axis. 
This  correction  is  illustrated  in  Figure  VI,   where  the  reader  should 
compare  the  curves  marked    d Q/  u  and    oO/L  (corrected),    but  has  been 
omitted  from  the  presentation  elsewhere  to  avoid  congesting  the  plots  un- 
necessarily.    In  the  first  cycle,   the  corrected  curve  gives    g**  '/  L  *gQ 

(At  this  point  in  all  hut  the  first  cycle,   each  value  of  shear  receives 
an  additive  correction  to  account  for  rotatory  inertia,   the  correction 
being  a  function  of  tae  bending  slope   7r  obtained  from  the  previous  cycle 

in  a  manner  to  be  described.     This  slope,   in  the  form  UY  ,   is  multiplied 

1 
by  the  working  function    ma  1,  L  A  and  the  product  added  to  the  shear 

term  to  give  a  curve  of 

tk&'/L,       m     gQ/u.     ♦      7|    tn    1  /   LA  (60) 
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This  is  the  only  procedural  respect  in  which  subsequent  cycles  differ 
from  the  first). 

The  curve  of    g  M'/  L  is  integrated,   and  the  uouadary  condition 
of  aero  end  moment  is  applied  ia  exactly  the  manner  described  auove 
for  the  shear  correction.     In  both  cases  the  correction  can  oe  considered 

to  filter  out  accumulated  errors  resulting  from  snail  inaccuracies  of 

2 

integraph  operation  and  calculation.     The  values  of     d  J*i   L   •   a  direct 

function  of  moment  in  ihe  hull  girder,   are  tabulated  from  the  corrected 

6    2 
curve,   and  each  is  multiplied  by  the  corresponding  value  of  10   L  /fcl. 

,      6 
to  give  a  function  of  curvature,        oV xlQ  ,   which  is  plotted. 

curvature  is  integrated  to  give  bending  slope  to  within  an 

additive  constant.     To  this  component  is  added  the  shear  component  of 

slope, 

6  6 

%/3x*j      L     «    (gQ/~)-r    (KAO/10   )  (61) 

based  on  the  shear  values  obtained  after  correction  of  the  first  integral 
in  the  cycle . 

The  addition  gives  a  function  of  total  slope,    still  without  the 
constant  of  integration: 

g7'  <  10   /L     u        <g*   *    i0   ,  ♦        g/3*  (62) 

To  determine  the  constant,   or  axis  shift,   use  is  made  of  Murray's 

formulation  (see  Part  I  of  this  appendix)  on  an  auxiliary  plot  (XXXV). 

6 
The  values  of       oyj'  «  10  .  L  from  the  above  addition  are  multiplied  by  the 

working  function 

and  the  product  plotted  and  integrated  on  the  auxiliary  plot.     The  value 
of  the  definite  integral  obtained, 
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is  added  to  the  uncorrected  total  slope  to  give  the  function  of  true  elope: 

6  6  6 

gy'xlO/i,      *        gy'xlO/^     +      lb  J  C./L  (64) 

This  is  plotted  on  the  cycle  plot  and  integrated,   yielding,   to 

within  an  additive  constant  of  integration,   a  function  of  deflection, 

6      2 

g  vf     x     10  ,  L   .     to  obtain  the  constant,   this  is  multiplied  by  mg 

and  the  product  plotted  on  the  auxiliary  plot    (XXIV)  and  integrated. 
The  value  of  the  definite  integral  is  divided  by  a  function  of  displace- 


ment to  give  the  constant 


0      4  A/  L 


6      2 
This  is  added  to  the  values  of    qyj  a  10      L     to  give    the 

deflection  pattern 

2  6      2  6      2 

OyxlO   /L        *   grjxlti  /l*        +    oC      *   10   i    . 

This  pattern  is  orthogonaiised  in  the  same  manner  as  was  des- 

6      2 
cribed  for  the  Kayieigh  pattern,   to  give     gy*  If  /L  ,     Ihc  value  at  the 

*       6      2 

normalising  point,        gy   «  10   /L   ,   is  recorded,   and  divided  into  all 

ordinate  values,   to  yield  the  normalised  deflection  pattern,        y  9 

of  the  first  cycle.     It  is  this  pattern  which  is  used  to  compute  the  inertia 

loading  for  the  next  cycle  of  the  iteration. 

6      2  2 

The  value  of      oy*      x     10  /L      is  multiplied  by  L   ,  and  the 

product  resulting  is  divided  into    the  values  of    fi&'x  10      from  the  last 
cycle  of  the  iteration  for  each  mode  to  yield  the  normalised  curvature 
for  the  mode,   which  now  will  have  the  correct  relationship  to  the 
normalised  deflection  pattern. 

Ihe  value  of   CO       corresponding  to  the  cycle  just  completed  is 

CO2      a L- —  (67) 
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2 

Ihe  value  of  cO       computed  from  the  last  cycle  performed  it  accepted 

as  the  natural  frequency  corresponding  to  the  mode  being  considered. 

In  preparation  for  the  rotatory  inertia  correction,   the  shear 

curvature  is  subtracted  from  the  total  curvature  to  give  the  corrected 

sending  curvature 

©6b 
gTftlQ   t  L     ■        g^'   x    10  ;  L    -      0(B     x      Iti  (feo) 

6      2 
which  is  normalised  by  dividing  by    oy    *     x    10   ft,      to  yield  the 

required  value  oi     ty  i,. 
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APPENDIX       t        -    SUMMAKY    OF    DMA    AND 


C  A  L  C  U  L  /  1  IONS 


Figuref    I  through    XXXIV    in  this  section  depict 
the  application  of  the  rr  ethod  developed  in  this  thesis 
to  the  model  for  which  characteristics  are  tabulated 
in  Appendix  t.     The  calculation  of  argument  curves 
plotted  from  data  is  illustrated  in  the  sample  calcula- 
tions oi  Appendix  C. 

Table  HI  summarizes  the  effect  of  each  iterative 
cycle,   and  its  associated  orthogonali&ation  check,   on 
the  natural  frequency  of  the  mode  considered.     The 
corresponding  effect  on  deflection  patterns  is  presented 
in  Pigu.rca  XII,   XX,   and    XXXII,   for  the  first,    second, 
and  third  modes  respectively. 
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Sample  Calculation 

'I his  appendix  contains  a  tabulation  of  the  working  .unction*  for 
iteration,   a  table  of  constants  used  in  the  iteration*   and  a  sample 
calculation  in  tabular  form.     All  symbols  arc  in  accordance  with  those 
appearing  in  the  Nomenclature. 

i  he  working  functions  (Table  IV)  pertain  to  all  oi  the  modal  calcu- 
lations.    These  functions-  were  calculated  as  a  matter  of  convenience  as 
well  as  necessity.     Their  form  is  governed  both  by  theory  and  t>y  the 
conventions  adopted  and  explained  in  i  ppendix  . 

I  able  V  contains  a  list  of  constants  which,   in  general,   apply  to 
all  three  modal  calculations.     Orthogonality  adjust  nent  constants  W    and 
W     apply  to  all  modal  calculations  beyond  the  first  mode  and  the  second 
mode  respectively.     1  he  derivation  and  explanation  of  the  constants  arising 
from  integration  may  be    ound  in  rtr:  1  of  appendix    >,     The  usual  value 
ot  Young's  Modulus  was  adopted,   and  the  shear  modulus  was  then  computed 
from  the  well  known  relation: 

c  ■  -iirr^r  • 

where  u    is  Pois son's  f\atio,    taken  equal  to  0.  3. 

The  sample  calculation  embraces  the  third  cycle  tl  the  third  mode. 
This  mode  was  chosen  a£     eing  the  most  general  of  the  three  with  respect 
to  the  orthogonality  conditions.     J. he  third  cycle  oi  this  mode  was  chosen 
tor  the  calculations  because  it  includes  the  rotatory  inertia  correction. 
(ihc  second  cycle  does  also,   but  the  results  thereo;  are  less  conclusive 
*han  those  of  the  third  cycle).     The  calculation  is  presented  in  tabluar  for a 
in   ..able  VI.     I  he  source  of  each  function  appears  beneath  its  symool  at 
the  head  of  the  column,     -o.^liiyiug  remarks  pertinent  to  the  sample 
calculation  appear  at  the  end  of  labie  VI. 
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'IAjbLJL    VI    (Continued) 

AMPCI*  YUNG  ftSMAMJ  1SLATIVI  rO  IAMPL1  CALCULATION 

*-- 

COLUMN  r-h-lmrtfcKS 

3    ..........   inertia  loading  with  CJ   taken  equal  to  1. 

5    ....... X.2.*1  tke  normalised  component  of  slope  due  to 

bending*   taken  from  the  previous  cycle. 

5,6 ?J.2  %"?£.     is  the  derivative  with  respect  to  x  of 

the  moment  of  the  inertia  forces,   at  any  section, 

about  the  vertical  centroidal  axis.     When  "K.t  o 

is  added  to  the  shear  force,  one  obtains     dM/dx, 
which  reflects  the  effect  of  lateral  inertia    orce* 
plus  shear  due  to  rotatory  inertia  within  the 
approximation  imposed  by  using  IT  from  the  previous 
cycle. 

9  ... ..... Component  of  slope  due  to  bending  moment  seiore 

determination  of  its  true  axis. 

10 . . Component  oi  slope  due  to  shear  flexibility. 

11  --.....-.... —  Column  9  plus  column  10.     Slope  before  true  axis 

determination. 
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17, 13, 19,  20  Calculation*  necessary  for  determination  of 

adjustment  constants  to  be  applied  to  — e~**3 
in  order  that  the  resultant  deflection  pattern  is 
orthogonal  with  respect  to: 

1.  igid  body  translation 

2.  Eigid  body  rotation 

3.  First  mode  deflection  pattern 

4.  Second  mode  deflection 

pattern 

17  1.    Orthogonality  with  respect  to  rigid  uody  translation: 

Condition:      /y3mdx    *    0 

csult: 
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lAfcJLi.  VI       Continued 


COLUMN 


18  2.     Orthogonality  with  respect  to  rigid  body  rotation 


Condition:    f%mU  (x  -  x)  dx    *    0 
•«> 

■eault: 

tm  « -^ — B l±JL 

-44.8    «    +0.1493 


300 
19  3.     Orthogonality  with  respect  to  Ut  mode: 


y6 


Condition:      /  y  .  y,  ro    dx    *    0 
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20  4.     Orthogonality  with  respect  to  2nd  mode 

Condition:         f y       y     m    dx     «    0 
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E  VI    Continued 


COLUMN  RI.MABKS 
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/*    6  2    - 

1    /  10  j    y    y 
/  "IT"  **  ' 

//*yy    y 


20  esult:      /  10  g  y    m    dx        a      -112 

m    dx 


w2 


-112  *     +    0.0222 


5040 


6 
d.iustroent  to  ue  applied  to    10    g    y  Cor 

L" 

orthogonality  with  first  mode  deflection  pattern. 


6      „ 

22  tauntf  to  Le  applied  to     10    g    y  for 

orthogonality  with  second    node  deflection,  pattern 

6     — 

23  Adjustment  to  i>e  applied  to     10    g    y  for 

orthogonality  with  rigid  body  rotation. 

24  J±My„*   iiJi-1,.3*    3  +  *<*-*>  ♦  Vi  +  CiiY2 

The  value  oi  10  g    ym  m  at  the  positive  anti-node 
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v       *** 


(near    X/l      »    0.5) 


is  the  normalizing  value  and  also  a  function  oi  the 
eigenvalue.     Its  value  If  J.  1503  for  this  cycle,   i.e., 
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c,  VI    Continued 


UMN  KS 


25  HwailiMl  resultant  deflection  pattern,  y 

>3»  ^  (0.1503) 


26  <on  of*/*   for  u#e  in  next  cycle. 
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APPI.Ni}U       D      -      ll»#JUgM»MTA»l    DISCUSSION 

1.     Secondary  Lffccts  Not  Treated: 

Shear  Cag;    Ihe  classical  beam  vU  rations  at  applied  to  ahipa 
do  not  take  special  notice  of  the  fact  that  resistance  to  deflection  is 
concentrated  generally  at  the  extremes  of  the  cross  section  and  that, 
in  horizontal  tending,   stress  transmitted  from  the  deck  to  the  side  will 
not  maintain  a  constant  level  down  the  side  as  indicated  ay  the  constant 
distance  from  the  n«utrai  axis.     Ihis  is  because  of  the  effect  known  as 
shear  lag.     While  Important  in  lightly  framed  ships  where  decks  and 
longitudinal  culkheads  play  little  part  structurally,   it  was  not  thought 
to  be  significant  in  the  ring-framed  model  of  this  investigation.     Shear 
lag  is  treated  in  detail  by  h.   A.   Anderson  and  J.   C.  Houoolt  ("The 
iLffect  of  Shear  Ltg  in  the  i.  ending  Vibrations  of  Box  Beams",   NACA 
TN  No.   1583,   May  1948)  and  others. 

Second  Frequency  Iff  tftl H    The  simple  Harmonic  assumption 
applied  to  iimoshenko's  complete  mechanism  leads  to  an  equation  which 
is  fourth  order  in  circular  iequency    ,     indicating  the  possibility  of  a 
second  real  value  of  frequency  for  a  given  :uode  .     Traill -Nash  and 
Collar     £63J     were  able  to  determine  this  second  Spectrum  of  frequencies 
for  certain  special  cases  of  uniform  beams.     Its  physical  importance  was 
not  determined,    out  it  •  was  thought  to  arise  from  a  resonant  interaction 
between  rotatory  inertia  forces  and  shear  stiffness. 

The  iterative  method  of  this  paper,  while  considering  both 
effects,  and  permitting  them  to  interact  in  the  sequence,  is  based  on  a 
single  assumed  frequency,  and  does  not  facilitate  simultaneous  inter- 
action.    Hence  the  second  frequency  spectrum  is  not  explicitly  noted. 

Torsional  Coupling:    It  should  be  noted  that  coupling  can  occur 
between  torsional  and  horizontal  flexural  vioratiou,  with  resultant  effects 
on  the  patterns  and  frequencies  of  each.     Such  coupling  has  not  I  sen  treated 
in  this  investigation. 
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Z.     More  Exact      rthogonality  Formulation: 

A  review  of  the  theoretical  development  of  this  method  revealed 
that  the  lornulatior.  of  the  orthogonality  relationships  used  were  inexact 
when  applied  to  orthogonali&ing  with  non-trivial  modes  when  rotatory 
inertia  is  considered.     The  effect  of  the  discrepancy  is  discussed  in  the 
next  su.  section.        The  development  of  the  more  exact  formulation* 
will  now  be  given. 

The  normalised  deflection  and  tending  slope,  and  the  circular 
frequency  of  each  of  two  different    nodes  are  given  -y 

yi*  *V    ^  I*       Yi    '    *J  *     ^J"  ' 
Multiplying  equations  (25)  and  (27)  I  endix  A-l    >y     If  -      and    V- 

respectively,   and  integrating,   we  have  for  mode   t    , 

|*j(M   ^')'dx       +     Ur-4yi*jdx     ♦        l«A«(y£/yj   -    Yt  ij  )dx     «    0   (69a) 

m  yt.  y.  co-  dx  +    /l**fl  ty.'  ■  *L )]  y-    •■    *   °  <**D> 

*  o 

and  using  the  multipliers  from  mode    i  ,  we  have  for  mode  j  , 

Ly./.w.dx    +     /Jkag(  y.    -  y.  ) 1  y.  dx       »    q  (70 .) 

Integrating  the  first  term  of  (69a)  ay  parts,   we  have 

.1  *•.*.'  ]     -  Jm  &}  dx    ♦  JlJ2~a>]*£*Jdx  t  /Loty/x-yx) 

0  dx      «    0  *  (71) 

♦Due  to  Prof.   5.  H.   Crandall,    Massachusetts  Institute  of  Technology 
and  to  C«    L»,   Dolph,    n  iilow  I«  research  Center,    University  oi  Michigan 
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where  the  first  term  vanishes  since  curvature  in  a  free-free  jeem 
vanishes  at  both  ends.     The  corresponding  operation  is  carried  out  on 
(70a),   and  the  result  subtracted  from  (71),   giving 


,^>|-te.r.7.«.  •fn-Wtj    -  1/Ti] 


dx      «     Q      (72) 


Now  (69b)  is  rewritten,  integrating  the  second  term  6     t  arts,   giving 

fm^6>?4*     ♦     KAG    ($'  *  *i  Wj]  ^^(y.'Vj.'  ^y/)dx-0    (73) 

where  the  second  term  vanishes,   since    KAG  (  y('  -    7^  )  is  identical  to 
shear,  which  must  vanish  at  both  ends  of  a  free -free  beam.     The    corres- 
ponding operation  is  carried  out  for  (70b),  and  the  resulting  equation  is 
subtracted  from  (73),   giving 

(  o)^  -  Wj*J    /      my-  f   dx       +  J  KAG  [%  f/    -tj  %']  dx      «     0  (74) 

Adding  (72)  and  (74)  gives  the  orthogonality  relationship 

(a)-1-  Uf)j[i-l2~7£r-  ♦      my.  y.    )    dx       -      I  (75) 

which,   for      co£  &  cJ-     ,   leaves  the  useful  relationship 
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/j 


This  should  be  compered  with  the  relationship  used  in  the  investigation 

(      y4   y-  )<**     ■  (77) 

which  is  seen  to  be  in  error  by  the  integral  of  the  first  term  in  (76). 

The  importance  of  this  ter-ri  varies  from    rode  to  mode,  and 
must  be  evaluated  for    higher  modes  to  determine  the  necessity  of 
applying  the  exact  condition.     It  does  not  lend  itself  readily  to  amplication 
in  the  iterative  procedure  of  this  thesis,   since  it  provides  corrections 
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to  bending  slop*  and  deflection,  quantities  which  oecome  available  at 
different  etagee  of  the  process,   and  which  cannot  M  corrected  separa- 
tely.    An  acceptable*   though  somewhat  unwieldy  approach  would  ue  to 
apply  the  complete  check  at  the  deflection  level,   examining  the  slope 
correction  to  see  if  its  magnitude  warrants  its  use.     If  it  is  significant, 
the  tending  elope  is  corrected,  and  the  new  total  slope  integrated.     The 
new  deflection  curve  is  given  the  exact  orthogonality  check,   this  time 
assuming  that  the  slo,  e  needs  no  correction. 

In  the  case  of  this  thesis,   an  auxiliary  investigation  showed 
that  the  complete  check  changed  the  natural  frequency  for  the  third  modes 
(previously  corrected  I     tne  approximate  relationship)  by  0.  36 /a,   and 
that  the  effect  was  even  smaller  on  the  second  mode  frequency,     it 
appears  that  for  the  modes  of  inter  est,   the  added  inconvenience  of  the 
full  correction  may  have  but  little  Justification.     Where  applied,   the 
correction  is  stated  fey  the  following  equations: 

tr.  ,c .  r.  *  *pr,    .    city  u  (79) 

where  the  constant  for  each  prior  mode  is  defined 
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3.     Third  S\ ode  Frequency  : 

As  noted  in  the  Discussion  of  Results,   the  natural  frequency 
obtained  for  the  third  mode  (1391    cpm)  _>y  the  iterative  process  and 
orthogonality  filtering  is  considered  to  be  slightly  high.     Ihis  stction 
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..resents  estimates  of  the  correct  value  obtained  from  several  independent 
approaches. 

a.     The  cycle -t,y~cycie  behavior  of  frequencies  in  the  third  mode 
calculation  manifests  alternation  casting  some  dou  t  on  the  value  from  the 
third  cycle.     Consideration  of  this  tehavior  seems  to  indicate  ultimate 
convergence  to  a  value  somewhere  between  the  second  and  third  cycle 
results,  whereas  earlier  modes  showed  continually  falling  trends.     The 
change  in  oehavior  in  the  third  mode  appears  to  stern  from  a  starting  as- 
sumption which  was  in  error,   even  after  orthogonali&ation,   to  a  greater 
degree  tnan  were  those  of  the  other  u*odes.     This  may  be  due  in  part 
to  the  use  of  the  inexact    orthogonality  check,   as  well  as  to  any  of  the  sources 
of  error  listed  in  the  Discussion  of  Kesults,     We  conclude  from  this  approach 
that  the  third  mode  natural  frequency  should     e  in  the  neighborhood  of 
1385  to  UN  cps. 

u.     The  more  exact  orthogonality  check  formulated  in  equation 
(§0)  was  applied  to  the  result  of  the  last  iteration    of  the  third  mode,     die- 
regarding  the  effect  on  slope.     1  he  constants     ecame 
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The  net  effect  was  to  lower  the  natural  frequency  from  1391  to  V3o6  tpm. 
The  change  in  n  ode  shape  was  insignificant. 
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v .     I  he  empirical  ratio  of  natural  frequencies  ooserved  for  a 
number  of    ships,   I:    2:    3,   although  approximate*  and  not  verified  for 
submarines,  is  applied  for  interest.     It  results  in  no  significant  change 
in  the  second  mode,    out  depresses  the  third  mode  frequency  to  1231 

>.     This  is  clearly    ouch  lower  than  there  is  reason  to  expect,    cut 
tends  to  verify  the  suspicion  that  1391  is  perhnpm  high. 

d.     Table  VII  lists  frequencies  for  the  model  under  considera- 
tion,  obtained  theoretically.     Tabulated  with  the  values  obtained  i~or*naUy 
in  this  investigation  (MIT)  are  those  ootained  by  David  Taylor  Model 
Basin  (DTMi)    y  analog  computer  and  those  obtained  oy  the  authors  at 
Underwater  Explosion  Research  Division  (Ut^O),    uoth  seta  including 

shear  flexibility  but  not  rotatory  inertia. 

41  VII 
CALCULATED  MQDI.L  FB&QUENCU.S  (cpn) 

Mode  MIT  DT.v  Ui-l^D 

I  417  444  42$ 

2  aaa  im  i»? 

3  1391  1392 

Since  the  added  refinement  of  rotatory  inertia  should  tend  to  lower  the 
frequencies  of  this  thesis  as  compared  to  the  other  sets,  it  again  appears 
that  1391  cpm  is  too  high  for  the  third  mode,  and  that  perhaps  a  reasonable 
value  to  expect  would  i.e  132  J  cpm.     This  is  based  on  both  the  ratio  series 
for  the  DTMi    data  and  the  suspected  trend  of  equal  or  greater  reduction 
with  mode  au.ia-er  due  to  rotatory  inertia. 

It  shodd  be  pointed  out  that  neither  set  of  comparison  frequencies 
is  corroborated,  and  experimental  frequencies  have  not  yet  been  obtained 
for  the  model  under  consideration. 

e.     Considering  all  of  the  above  factors,   it  can  best  be  concluded 
that  the  third     .ode  frequency  is    1355    +  2.6  k. 
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APPENDIX         £         -         STRUCTURAL    DATA 

Table  VIII  contains  a  tabulation  of  structural  data  necessary  for 
the  modal    calculations.     The  source  of  these  data  is  the  Underwater 
Explosions  Research  Division    of  the  Norfolk  Naval  Shipyard,   under  whose 
auspices  the  model  submarine  referred  to  in  this  work  was  constructed. 

The  mass  appearing  in  Table  VIII  includes  the  structural  mass 
of  the  submarine,   the  mass  of  her  ballast  water,    the  mass  of  her  solid 
contents,   and  the  virtual  water  mass. 

The  virtual  water  mass  was  calculated  from  the  following 
formula: 

mv    =  £7T  b  (81) 

where    b    is  half  the  depth  of  the  submarine  structure  at  any  section, 
including  the  protuberance  of  the  rudder  and  conning  tower  where  appli- 
cable,   and     f.     density  of  water  in  lb/ft3. 

Equation  (81)  is  based  on  the  assumption  that  the  cross  section 
is  elliptical  with  the  major  axis  congruent  with  the  vertical  center  line 
of  the  submarine.     The  wash  of  the  water  through  the  superstructure  is 
neglected.     The  resulting  curve  of  virtual  water  mass  was  further  modi- 
fied by  an  empirical  correction  explained  in  David  Taylor  Model  Basin 
Report    #   632.     This  correction  consisted  of  multiplying  the    mv   curve 
by  0.  83  over  most  of  its  length,   with  a  more  severe  correction  at  the 
ends.     Only  one-third  of  the  additional  virtual  water  mass  contributed 
by  the  rudder  and  the  conning  tower  was  used. 

These  corrections,   am    somewhat  arbitrary,   as  all  estimates  of 
virtual  water  mass  tend  to  be.     Tne  total  empirical  correction  causes 
only  a  very  slight  change  in  the  ship  vibration  behavior.     No  experimen- 
tal data  that  verifies  the  virtual  water  mass  distribution  is  available 
for  this  submarine. 

The  ends  of  the  total  mass  curve  used  oy  the  authors  were  altered 
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•lightly  so  that  the  ordinates  did  not  approach  aero  rapidly,    out  remained 
at  an  ar  itrary  small  amount,   dropping  abruptly  to  zero  at   *  *    0     and 

X  ■      L,     This  was  done  to  conform  with  a  similar  adjustment  to  the  mom- 
ent of  inertia  curve.     The  need  for  the  corrections  arose  when,  in  their 
absence,    spurious  values  of  curvature  were  obtained  at  the  ends  because, 
in  general,   the  curve  of  inertia  tended  to  converge  on  zero  more  rapidly 
than  the  curve  of  bending  moment. 

The  K    factor  used  to  determine    O    from      Q     is  defined  as 

average  shearing  stress  in  cross  section 
maximum  shearing  stress  in  cross  section 

or 

V  

*      A    C\    y</A  (82) 


a  r\ 


where  I    and  A    have  already  been  defined,  and  B    is  the  total  thickness 
of  steel  through  which  the  vertical  center  line  passes.     (The    center  line 
vertical  and  flat  keels  were  excluded  from  the  calculation  of  K). 

The  integral  of  equation  ($2)was  calculated  separately  for  each 
component  of  the  structure  (i.e.,   the  pressure  hull,   outer  hull,   deck, 
etc.  ),  and  all  the  integrals  summed  to  give  a  value  applicable  to  the 
entire  section.     Only  members  which  make  rr.sjor  contributions  to  the 
submarine's  longitudinal   strength  were  considered  in  the  above  calcu- 
lation and  in  computing  sectional  moments  of  inertia.     A  discussion  as 
to  the  validity  of  this  practice  appears  in  Appendix  D. 
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